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QCD with dynamical Wilson fermions at /3 = 5.5* 

K.M. Bitar, R.G. Edwards, U.M. Heller, and A.D. Kennedy'^ 
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We study QCD with two flavors of dynamical Wilson fermions at /3 = 5.5 and three values of k. The corre- 
sponding pion masses are 0.375, 0.324 and 0.262 in lattice units, with pion to rho mass ratios of 0.76, 0.71 and 
0.62, respectively. We use the configurations to compute the heavy quark potential, leading to lattice spacings of 
0.110, 0.105 and 0.099 fm, and to compute spectroscopy for several different valence quark k's. 
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1. INTRODUCTION 

Of the two popular unimproved ways of dis- 
cretizing the Dirac operator simulations with dy- 
namical Wilson fermions still lag behind simu- 
lations with dynamical staggered fermions. The 
best completed simulation with Wilson fermions, 
by the HEMCGC collaboration |]l|, was done at 
stronger couplings, i.e., larger lattice spacing, 
than the state-of-the-art staggered simulations 
|2||^, and at larger m^/wp ratios, 0.72 and 0.60. 
The lowest m^/rTip ratio achieved to date with 
staggered fermions, in contrast, is 0.33 Q, albeit 
at quite a strong coupling. There are many indi- 
cations that the Wilson fermion HEMCGC sim- 
ulations were done far from the continuum limit: 
for example a~^, in the chiral limit, was estimated 
as 1600 MeV from the nucleon mass, 1800 MeV 
from the rho mass and 2000 MeV from , and a 
computation of A/3 proved impossible [|| . 

About 3 years ago, the SCRI lattice group de- 
cided to push dynamical Wilson fermion simula- 
tions to a somewhat weaker coupling with runs 
on the SCRI CM-2. Being limited, for technical 
reasons, to lattice sizes that are a power of two, 
we chose a 16^^ x 32 lattice and a gauge coupling 
P = 5.5. With these choices we hoped that the fi- 
nite size effects will be relatively small. We have 
made runs at three different k values, to allow 
chiral fits and thus to test the chiral extrapola- 
tions. In each case we have about 2000 trajecto- 
ries, after thermalization, with parameters given 
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Table 1 

Summary of run parameters 



K 


dt 




R 


Acc 


(Ncg) 


0.1596 


0.01754 


3 


10-' 


75% 


270 


0.1600 


0.01515 


3 


10-7 


75% 


380 


0.1604 


0.007 


3 


10-7 


90% 


450 



Table 2 
Some results 



K 






^plaq 


T7r(10) 


0.1596 


0.55928(9) 


0.178 


33 


65 


0.1600 


0.56025(5) 


0.177 


9 


15 


0.1604 


0.56131(8) 


0.176 


21 


21 



in Table Table ^ contains some results, includ- 
ing rough estimates of integrated autocorrelation 
times. Sample time histories, for the run with 
K = 0.1600, are shown in Fig. |l|. 

2. THE HEAVY QUARK POTENTIAL 

On the ensembles of configurations generated, 
we computed the heavy quark potential from a 
combination of on- and off-axis timelike Wilson 
loops, with "APE smeared" links in the space 
directions |^ to enhance the signal to noise ra- 
tio. Since string breaking is estimated to occur 
at a distance somewhat larger than accessible to 
us (and we do not observe any sign of it) we used 
the same ansatz that is common in quenched sim- 
ulations for fitting the potential 
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Figure 1. Time history of the plaquette and the 
pion correlation function at distance 10 with waU 
source at to = and 10 for k — 0.1600. 
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Figure 2. The potential, in physical units, includ- 
ing, for comparison, the (3 — 5.3 data from M. 



with Gl, the lattice Coulomb potential, taking 
account of the lattice artefacts present at small 
distances. We used fully correlated fits with 
the covariance matrix estimated by a jackknife 
method. 

The parameters for the best fits, defined as in 
ref. 0], are listed in Table ||. We show the po- 
tentials together with the fits, in physical units 
as determined from the string tension, in Fig. |^. 
Also included in the plot are the potentials from 
the (3 = 5.3 simulations as computed in ref. |^]. 
All curves agree fairly well, showing approximate 
scaling of the potential. The Coulomb coefficient 
e is somewhat larger than for quenched calcula- 
tions at comparable lattice spacing, which can be 
attributed to the effect of the dynamical fermions 

i- 

3. SPECTROSCOPY 

For each ensemble we computed the hadron 
spectroscopy with six different valence quark kap- 
pas, Kv = Ksca, 0.1590, 0.1580, 0.1395, 0.1230 and 
0.1050. The last three are quite heavy, in the c 
to b quark range, and will not be discussed here. 
Mesons were computed for all pairs of kv 's, while 
baryons only for equal valence quarks. Mass fits, 
preferably to two states, and chiral fits were fully 



correlated. In Fig. |3| we show chiral plots (vs. 
1/k) for the results with Kgoa = 0.1600. Similar 
plots for the other sea quarks and tables of all the 
masses can be found in ref. . 

Having ensembles with three different sea 
quark masses, we can also make chiral fits for 
spectroscopy involving only Ky = Ksqh- These are 
shown in Fig. ^. We obtain Kc — 0.16116(15) and 
extrapolated masses anip — 0.352(12), amj\j = 
0.520(18) and otoa = 0.603(15). These numbers, 
as well as the results in Table ||for Kgca = 0.1604, 
are still preliminary. We are in the process of 
computing propagators from another time slice. 

In Table ^ we give a collection of estimates of 
in MeV, obtained from the potential, rho, nu- 
cleon and Delta masses. For the latter, at fixed 
Kscn, the masses were extrapolated to the physi- 
cal m-jj /mp ratio. The last line shows the extrap- 
olation to zero sea quark mass from the masses 
obtained with Ky — Ksca- As can be seen, the 
lattice spacing estimates still differ considerably 
across each line in the table and we have to con- 
clude that there are still rather large 0(a) effects 
present. It will be difficult to decrease these lat- 
tice effects by simply going to smaller lattice spac- 
ings. 
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Table 3 

Sximmary of results from fits to the effective potentials using eq. (|l|) . The last two columns give the scale 
ro/a determined from rgi^(ro) = 1.65 and the dimensionless quantity r^y/a. 



K 


T 


^min ^max 


aVo a^<7 e e' 


Q 




0.1596 


4 


1.41 13.86 


0.731(3) 0.0645(7) 0.315(3) 0.368(12) 


0.90 


4.55(2) 1.155(1) 


0.1600 


4 


1.73 13.86 


0.743(5) 0.0583(9) 0.329(6) 0.375(13) 


0.77 


4.76(3) 1.149(3) 


0.1604 


4 


2.24 13.86 


0.757(6) 0.0510(9) 0.342(5) 0.53(8) 


0.81 


5.06(3) 1.144(4) 
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Figure 3. Chiral plots for the hadron spec- 
troscopy on the ensemble with Kgca = 0.1600. The 
dashed lines indicate Kc- 
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Figure 4. Chiral plots for the hadron spec- 
troscopy with Ky ~ Ksoa- 



Table 4 

Values for a^^ in MeV. The last line comes from 
extrapolations in the sea quark mass to the chiral 
limit. 





p N A 


0.1596 
0.1600 
0.1604 


1790(10) 1955(40) 1650(30) 1910(45) 
1885(15) 2025(25) 1685(30) 1950(45) 
1995(15) 1995(40) 1755(45) 1985(40) 


0.16116 


2185(25) 2190(75) 1805(65) 2045(55) 
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